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In this paper we analyze the impact on the orbital motions of the outer planets of the 
solar system from Jupiter to Pluto of some velocity-dependent forces recently proposed 
to phenomenologically explain the Pioneer anomaly, and compare their predictions (sec- 
ular variations of the longitude of perihelion xa or of the semimajor axis a and the 
eccentricity e) with the latest observational determinations by E.V. Pitjeva with the 
EPM2006 ephemerides. It turns out that while the predicted centennial shifts of a are 
so huge that they would have been easily detected for all planets with the exception of 
Neptune, the predicted anomalous precessions of w are too small, with the exception of 
Jupiter, so that they are still compatible with the estimated corrections to the standard 
Newton-Einstein perihelion precessions. As a consequence, we incline to discard those 
extra-forces predicting secular variations of a and e, also for some other reasons, and 
to give a chance, at least observationally, to those models predicting still undetectable 
perihelion precessions. Of course, adequate theoretical foundations for them should be 
found. 

Keywords: Experimental tests of gravitational theories; Modified theories of gravity; 
Celestial mechanics 



1. Introduction 

The Pioneer anomaly^ (PA) consists of an unmodelled almost constant and uniform 
acceleration approximately directed towards the Sun of magnitude 

A Pio = (8.74 ± 1.33) x lO" 10 m s" 2 (1) 

detected^™ in the radiometric data from the Pioneer 10 (launched in March 
1972) and Pioneer 11 (launched in April 1973) spacecraft after they passed the w 20 
AU threshold moving with speed vpi « 1.2 x 10 4 m s _1 along roughly antiparallcl 
escape hyperbolic paths undertaken after their previous encounters with Jupiter 
(w 5 AU) and Saturn (ps 10 AU), respectively. Concerni ng t he possibility that PA 
started to manifest itself at shorter heliocentric distances^! efforts t o re trieve and 
analyze early data from Pioneer 10/11 are currently being performecPM 
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The Pioneer spacecraft were particularly well suited for radioscience celestial 
mechanics experiments because they were spin-stabilized^; in practice, they could 
be regarded as gyroscopes so that only a few orientation maneuvers, easily modeled, 
were needed every year to keep the antenna pointed towards the Earth. On the 
contrary, 3-axis stabilized spacecraft like Voyager 1/2 undergo continuous, semi- 
autonomous, small gas jet thrusts to maintain the antenna facing the Earth; as a 
consequence, their navigation is not as precise as that of the Pioneer 10/11. 

The attempts performed so far to ex plain PA in terms of known effects 
of gravitat ional^ and/or non-gravitationa l^^ ^ origin were found to be not 
satisfactor y^ 2 | 13 |^ g0 ^ na ^- a vag £ numDer f exotic explanations based on modified 
models of gravity were proposed (see, e.g., Ref. 151 [T H \15 \ [TB I an d references therein). 
If PA is due to some modifications of the known laws of gravity, this must be due 
to a radial extra-force affecting the orbits of the planets as well, especially those 
moving in the region in which PA manifested itself in its presently known form. The 
impact of a Pioneer-like additional acceleration on the motion of major and minor 
bodies in the outer regions of th e solar system wa s recently studied by numerous 
authors with different a 

pproache Jl2llHI19l2Ql2Tj22l 

it turned out that a constant 
and uniform extra-acceleration with the magnitude of Eq. ([1]) would produce huge 
secular effects which are ne atly absent in the planetary data. 



It was recently suggested^ that, from a purely phenomenological point of view, 
test bodies moving in the (outer) solar system could experience velocity-dependent 
extra-accelerations of the form 



where v r is the radial component of the test particle's velocity v; Eq. (J2J) and Eq. 
Q would reduce to Eq. |T]) for the Pioneer 10/11 spacecraft whose velocities can 
be assumed entirely radial in the outer regions of the solar system in which PA 
was detected. Standish in Ref. put on the test such a hypothesis by fitting huge 
planetary data sets with the dynamical force models of the latest Jet Propulsion 
Laboratory (JPL) DE ephemerides modified ad hoc according to Eq. (J2J) and Eq. 
(|3"|) and examining the results in terms, e.g, of the reliability of the estimated pa- 
rameters. His conclusion was that the existence of extra-accelerations like those of 
Eq. @ and Eq. ^ at heliocentric distances > 20 AU cannot be ruled out by the 
present-day available data of the outer planets because Eq. @ and, especially, Eq. 
([3]) would induce orbital effects on them too small to be detected. Their existence 
in the inner regions of solar system is, instead, ruled out. 






and 





a This was due to the fact that Pioneer 10/11 were equipped with Radioisotope Thermoelectric 
Generators (RTG) placed at the end of long booms to be away from the spacecraft and thereby 
avoid any radiation damage. 
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In this paper we will follow a different approach by using the EPM2006 
ephemerides produced by E.V. PitjevaPSl at the Institute of Applied Astronomy 
(IAA) of the Russian Academy of Sciences (RAS). First, we will analytically work 
out the secular effects of small perturbing accelerations like those of Eq. ((2|) and 
Eq. ([3]) on the Keplerian orbital elements of a planet in order to gain as clear as 
possible insights about the modifications which the orbits would undergo if Eq. ((2|) 
and Eq. ((3]) were real; should some implausible physical feature turn out, it would 
be more difficult to trust such proposed anomalous forces. Then, we will compare 
some of such predictions with the latest observational determinations for the outer 
planets estimated by Pitjeva with the EPM2006 ephemerides in a purely phemo- 
menological way as corrections to the known effects due to usual Newton-Einstein 
laws, without modelling any additional force. In Table Q] we quote some quantities 
we will use. They are the outco me of a global fit of more than 400,000 data points 
(1913-2006) performed by PitjevgPSEl with the EPM2006 ephemerides; about 230 
parameters were estimated. It must be noted that the uncertainties 5 Am in the 

Table 1. Second column: formal standard deviations <5a, in m, 
of the semimajor axes of the outer planets from a fit of 400,000 
data points spanning almost one century with the EPM2006 
ephemerides (from Table 3 of Ref. 24). Third column: correc- 
tions to the standard Newton-Einstein secular precessions of the 
perihelia 25 in arcsec cy — 1 . Fourth column: their formal errors, 
in arcsec cy — , re-scaled by a factor 10. For Neptune and Pluto 
no secular precessions have been estimated because the available 
modern data records for them do not yet cover an entire orbital 



revolution. 


Planet 


8a (m) 


A to (arcsec cy 1 ) 


<5Aro (arcsec cy 1 ) 


Jupiter 


615 


0.0062 


0.036 


Saturn 


4,256 


-0.92 


2.9 


Uranus 


40,294 


0.57 


13.0 


Neptune 


463,307 


N.A. 


N.A. 


Pluto 


3,412,734 


N.A. 


N.A. 



estimated corrections to the perihelion precessions are the formal ones re-scaled by 
a factor 10 in order to obtain realistic evaluations for them. 

2. The orbital effects of velocity-dependent perturbing forces 
yielding Pioneer-type accelerations 

2.1. Forces linear in velocity 

According to the classification of Ref. [22j the first two kinds of extra-forces are 
linear in v r being 

4 2) = -k|/c, (4) 

4 3) = -VrK, (5) 
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with 

K = — — = 7.3 x 10" 14 rad s" 1 = 47.4 arcsec cy" 1 . (6) 

Wpio 

The radial acceleration of Eq. is constantly inward, i.e. directed towards the 
Sun, while the one of Eq. (f5|) is directed towards the Sun when v r > 0, i.e. when 
the planets gets farther from the Sun, while is directed away the Sun when v r < 0, 
i.e. when the planet gets closer to the Sun. Indeed, for an unperturbed Keplerian 
ellipse^] 

nae sin / nae sin E 



y/1 - e 2 1-ecosS' 



(7) 



where a is the semimajor axis, e is the eccentricity, n = yj GMj a 3 is the mean 
motion, / is the true anomaly counted anticlockwise from the perihelion, and E 
is the eccentric anomaly. v r > for < / < 7r, i.e., from the perihelion to the 
aphelion, and v r < for tt < f < 2ir, i.e. from the aphelion back to the perihelion. 

In view of the smallness of Eq. (j4|) and Eq. ([5]) for the planets of the solar system, 
we will treat them perturbatively. Indeed, the radial velocities for the outer planets 
amount to 10 1 — 10 3 m s _1 only, so that A v as 10 -11 m s -2 , while the Newtonian 
attraction of the Sun is for them of the order of 10~ 4 — 10~ 6 m s~ 2 . 

Let us work out the secular precession of the longitude of perihelion w. The 
Gauss equation for its variation due to an entirely radial perturbing acceleration 
A r id™ 



dw y/1 - e 2 

— = A- cos/. (8 

at nae 

By inserting Eq. ([4]) into Eq. |(5J), evaluating the r.h.s over the unperturbed Keplerian 
ellipse characterized by 

r = a(l — e cosE), 
dt = ( 1 ~ e T '; os - E ) dE, 
cos / = -f-° sE ~ 



(9) 



-c cos 



E' 



sin/= ' / F^ S % E , 

J 1-ecosE 1 

and averaging over one orbital period one gets 

^7-73 ' 2 ^t^ 

(ZD) = 



7T 



2e-(l-e 2 )lnU±f 



<0. (10) 



We have used the fact that 



\v r \ = v r , < / < 7r; \v r \ = — v r , 7r < / < 2tt. 



(11) 
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Instead, Eq. (JSJ) yields no perihelion precession. Indeed 
ICVT~^ f 2 * (cosE~e)smE 



dE = 0. (12) 
2tt J l-ecos£ v ; 

The Gauss equation for t he variation of the semimajor axis due to the a radial 



perturbing acceleration is^l 

da 2e . , . 

^ = WT=^^ sin/ - (13) 

By proceeding as before it turns out that Eq. (H|) does not yield secular variations 
of a; instead, Eq. (|5|) induces a secular decrease of a according to 

(d) = 2/Ca (^=L= - l) . (14) 

Note that for circular orbits, i.e. v r = 0, (d) = 0. Since 

de _ ( 1 — e 2 \ da 
dt ~ \ 2ae J ~dl 

when A = A r , also the eccentricity decreases: 



(15) 



1 -1). (16) 



As expected for a central force, the orbital angular momentum L = V GMa(l — e 2 ) 
is conserved, on average: indeed, Eq. (fT4"|) and Eq. (fll)|) yield 

(^)-GAf[(d)(l-e 2 )-2 fle (e)] =0. (17) 

Instead, the energy £ = —GM/2a is not conserved; according to Eq. (fT4| . 
• \ GM ... GMK (I \ n 

Such a result is certainly suspect from a physical point of view. 

In order to independently check the results obtained analytically we performed 
two numerical integrations of the equations of motion adding to the Newtonian 
monopole term the perturbing accelerations of Eq. (j4j and Eq. ((5]) : the qualitative 
features of the resulting motions are depicted in Figure [T] and Figure [2] 

Let us now consider the problem of the existence of the accelerations of Eq. (j4|) 
and Eq. ([5|) from a phenomenological point of view according to the present-day 
planetary data available. In Table [2] we quote the predictions for the outer planets 
of the centennial shifts in m of the semimajor axis, according to Eq. (fT4|) . and of the 
secular perihelion precessions in arcsec cy _1 , according to Eq. (flQ|) . Such predictions 
must be compared with the observationally determined parameters quoted in Table 
[TJ Concerning the semimajor axis, the present-day accuracy in determining them 
would clearly allow to detect shifts as large as those of Table [2] for all planets from 
Jupiter to Pluto with the exception of Neptune, even by re-scaling the values of 
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y, AU 




Fig. 1. Numerically integrated trajectory for the radial acceleration of Eq. J4} proportional to 
— \v r | . The longitude of perihelion zu undergoes a retrograde precession while neither the semimaj or 
axis a nor the eccentricity e experience secular variations. 

y, AU 




Fig. 2. Numerically integrated trajectory for the radial acceleration of Eq. ([5} proportional to 
— IV- Both the semimajor axis a and the eccentricity e secularly decrease, while the longitude of 
perihelion xa remains fixed. 

Table 2. Second column: shift Aa, in m, 
of the semimajor axis of the outer plan- 
ets over 1 cy according to Eq. (1146 induced 
by the acceleration of Eq. (0 proportional 
to — v r . Third column: secular precessions 
■dj of the perihelia of the outer planets, 
in arcsec cy , according to Eq. I|10[ l due 
to the acceleration proportional to — |iv|. 
Such values are to be compared with those 



in Table Q] 


Planet 


Aa (m) 


w (arcsec cy 1 ) 


Jupiter 


-419,726 


-0.973 


Saturn 


-1,030,797 


-1.1 


Uranus 


-1,470,544 


-0.9 


Neptune 


-76,220 


-0.1 


Pluto 


-88,154,057 


-4.9 



5a of Table Q] by a factor 10 or more. The situation is less neat for the perihelion 
precessions. Indeed, it turns out that the present-day accuracy in determining them 
does not allow to rule out Eq. (jlOjl . with exception of Jupiter. Thus, we conclude 
that the acceleration of Eq. (HJ proportional to — v r is to be considered ruled out 
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by observations; it is true that, in principle, adjusting the ephemerides (without 
modifying their dynamical force models as done with EPM2006) may absorb the 
exotic signatures, but we do not believe that could occur because of the huge size 
of them. Instead, the effects induced by Eq. (|4]), proportional to — \v r \, are still 
compatible with data. Of course, in drawing such conclusions we are tacitly assuming 
that the Pioneer-type anomalous accelerations of Eq. (|4]) and Eq. (O exist since one 
century at least. Let us assume that they act since much longer time, say 500 Myr; 
Eq. ([5]) and Eq. (fT4"]) tell us that, in this case, 500 Myr ago the semimajor axes of 
the outer planets were equal to 19 AU (Jupiter), 44 AU (Saturn), 68 AU (Uranus), 
32 AU (Neptune) and 2,983 AU (Pluto). We have used the simple formula 

ao = a — aAt, (19) 

in which ao represents the semimajor axis in the past while a denotes its current 
value. Concerning the eccentricities, they would have been larger than 1 according 
to Eq. (|T6|) and eo = e — eAt. With regard to the future evolution of the orbits of 
the outer planets, the time required to circularize their orbits with respect to the 
present-day values of the eccentricities is of the order of 8 x 10 5 yr, provided that 
the Pioneer-type forces considered here will continuously act upon the planets for a 
so long time span. Of course, issues concerning a theoretical justification for Eq. (j4|) 
and Eq. (|5|) remain: suffices it to say that they are not, in general, Lorentz-invariant, 
as can be straightforwardly shown by using 

r = F(r - Vt) + (T - l) ^^ , t = r ( t - , (20) 



1 vV 





/r-n 


v- V + 




V TV 2 J 



(v x V) x V 



(21) 



with T = - V 2 /c 2 . The conclusions by StandisrP^ 1 are that Eq. |@J and 

Eq. ([5]) cannot exist for planets up to Jupiter and Saturn, while their existence at 
heliocentric distances > 20 AU is virtually undetectable from the motion of Uranus, 
Neptune and Pluto. 



2.2. Forces quadratic in velocity 

The other two anomalous radial accelerations examined in Ref. 1221 quadratic in the 
radial velocity, are 

A$ = -v 2 ,H, (22) 
A$ = -\v r \v r H, (23) 

with 

W = ^ = 6 .07xl0- 18 m- 1 . (24) 

V Pio 
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The acceleration of Eq. (|22p is always directed towards the Sun, while the one of 
Eq. (|23|) is inward when the planet moves away from the Sun, while it is directed 
outwards when the planet approaches the Sun, as in the case of Eq. ([5]) . 

An acceleration like Eq. (|22l) was theoretically obtained by Jaekel and Reynaud 



in the framework of their lineaJ^-^H] anc j non-lineaJ22l metric extensions of general 
relativity. Its orbital effects were worked out in Ref. [T9j neither the semimajor axis 
nor the eccentricity undergo secular variations, while the longitude of perihelion 
precesses according to 



HnavL 



2\/l 



< 0. 



(25) 



In Figure [3] we show the results of the numerical integration of the equations of 
motion with Eq. ([22| added to the Newtonian monopole: the results obtained ana- 
lytically in Ref. QjJ] are confirmed. 




Fig. 3. Numerically integrated trajectory for the radial acceleration of Eq. 1221 ) proportional to 
— v%. Both the semimajor axis a and the eccentricity e remain unchanged, while the longitude of 
perihelion zu experiences a retrograde secular precession. 



In Ref. 1191 it was shown that the inner planets' perihelion precessions predicted 
by Eq. (|25[) are neatly ruled out by the corrections to the perihelion precessions 
estimated by Pitjeva^ with the EPM2004 ephemerides. In Table [3] we quote the 
predictions for the outer planets; it turns out that they are compatible with the 
results of Table [IJ apart from Jupiter. Note that it is true also by considering the 
formal uncertainties in the estimated corrections to the perihelion precessions, i.e. 
the values of 5Avj in Table [1] re duced by 10 times. Such a conclusion substantially 
agrees with that by Standisl ^Zl. 

Eq. (|23|) . contrary to Eq. I|22[). induces no secular perihelion precession and a 
secular variation of the semimajor axis and the eccentricity which decrease according 
to 



4Hna 2 



2e + ln 



(26) 
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Table 3. Secular preces- 
sions w of the perihelia of 
the outer planets, in arc- 
sec cy — , according to Eq. 
(125 I t due to the acceleration 
of Eq. \22\ proportional to 
— v^. Such values are to be 
compared with those in Ta- 
ble [T] 



Planet 


uj (arcsec cy 1 ) 


Jupiter 


-0.030 


Saturn 


-0.029 


Uranus 


-0.015 


Neptune 


-0.0004 


Pluto 


-0.308 



(£)]■ 

Note that (a) — for circular orbits. Thus, also Eq. (|23|) does not conserve the 
total energy. Such analytical results are confirmed by a numerical integration of the 
equations of motion showed in Figure [4j 



2Una{\ - e z ) 




Fig. 4. Numerically integrated trajectory for the radial acceleration of Eq. (12:{|i proportional to 
— |tv|?V. Both the semimajor axis a and the eccentricity e secularly decrease, while the longitude 
of perihelion zu remains unchanged. 



In Table[4]we quote the predictions for the centennial semimajor shifts according 
to Eq. (|26| . A comparison with the Table [1] shows that the formal uncertainties 8a 
are always quite smaller than such anomalous shifts, apart from Neptune. However, 
it must taken into account that realistic errors may be up to one order of magnitude 
larger: if so, it would not be possible to rule out the results of Table [4] apart 
from Jupiter. In this case, our conclusions would agree with those by StandisrP^J 
Of course, serious issues concerning theoretical justifications of Eq. ([23| and the 
temporal extent of its existence remain open. Indeed, given the present-day values of 
the planetary semimajor axes and eccentricities and assuming that Eq. ()23|) existed 
unchanged in the deep past, about 100 Myr-lGyr ago e = 1 for the planets from 
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Table 4. Shifts of 
the semimajor axes 
a of the outer plan- 
ets, in m, according 
to the acceleration of 
Eq. I|23j l proportional 
to —|n r |iv. Such val- 
ues are to be com- 
pared with those in 
Tabled] 



Planet 


Aa (m) 


Jupiter 


-18,753 


Saturn 


-39,362 


Uranus 


-33,347 


Neptune 


-251 


Pluto 


-7,283,499 



Jupiter to Pluto. Since, instead, the semimajor axes would have remained almost 
unchanged, this means that the perihelion distances vanished. 

3. Conclusions 

An ingenious attempt recently proposed to explain the Pioneer anomaly as due to a 
modification of the usual Newton-Einstein laws of gravitation consists in postulat- 
ing the existence of some velocity-dependent extra-forces linear or quadratic in the 
radial component v r of the velocity of a test body. We put on the test such empirical 
models in the outer regions of the solar system in which the Pioneer anomaly mani- 
fested itself in its presently known form with the latest observational determinations 
of the planetary motions obtained by E.V. Pitjeva with the EPM2006 ephemerides. 
It turns out that the models yielding anomalous perihelion precessions cannot yet 
be ruled out, at least phenomenologically, for heliocentric distances larger than 5 
AU. On the contrary, the models predicting secular variations of the semimajor 
axis a and the eccentricity e are much more difficult to be trusted not only because 
they would violate the conservation of energy but also because the centennial shifts 
for a predicted by them are so large that they should have been detected, given 
the present-day accuracy in determining such orbital element. However, it must be 
considered that sound theoretical justifications for such models must be given. 
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